Mass dispersion in transfer reactions with a stochastic mean-field theory 
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Nucleon transfer in symmetric heavy-ion reactions at energies below the Coulomb barrier is in- 
vestigated in the framework of a microscopic stochastic mean-field theory. While mean-field alone is 
known to significantly underpredict the dispersion of the fragment mass distribution, a considerable 
enhancement of the dispersion is obtained in the stochastic mean-field theory. The variance of the 
fragment mass distribution deduced from the stochastic theory scales with the number of exchanged 
nucleon. Therefore, the new approach provides the first fully microscopic theory consistent with the 
phenomenological analysis of the experimental data. 
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The mean-field theory, otherwise known as the time- 
dependent Hartree-Fock (TDHF) model [ll, i, i, i, i, 0] , 
has been widely used in describing reaction dynamics at 
low energies in nuclear collisions and other many-body 
systems. It is well known that the mean-field theory pro- 
vides a good description of mean values of one-body ob- 
servables in low energy reactions. However, it completely 
fails in the description of the dynamics of fluctuations of 
one-body observables. During the past decades, large ef- 
forts have been devoted to overcome this difficulty and to 
develop transport theories that are able to describe not 
only mean values but also fiuctuations (for a review see 
0, [§])■ Among them, the variational principle proposed 
by Balian and Veneroni (BV) appears as one of the most 
promising methods [ol [lol . [ill . However, even nowadays 
it remains difficult to applv [l3. ITsj . More than 30 years 
after the first application of TDHF, the absence of prac- 
tical solution to include fluctuations beyond mean-fleld 
in a fully microscopic framework strongly restricts appli- 
cations of mean-field based theories. 

There arc mainly two mechanisms for density fluctua- 
tions: (i) collisional fluctuations generated by two-body 
collisions and (ii) one-body mechanism or mean-fleld fluc- 
tuations. At low bombarding energies, the mean-fleld 
fiuctuations provide the dominant mechanism for fiuctu- 
ations of collective motion and collisional effects could be 
neglected. Restricting the treatment at low energies, re- 
cently, a stochastic mean-field (SMF) approach has been 
proposed for nuclear dynamics [l3| . It was demonstrated 
that the approach incorporates the one-body dissipa- 
tion and associated fiuctuation mechanism in accordance 
with the quantal dissipation-fiuctuation relation. Fur- 
thermore, in the limit of small amplitude fluctuations, 
the SMF approach gives the same result for dispersion 
of one-body observables as that of the formula derived 
from the BV approach JLQJ. Therefore, the SMF ap- 
proach provides a powerful tool for describing low energy 
nuclear processes including induced flssion, heavy-ion fu- 
sion near barrier energies, dccp-inclastic collisions, and 
spinodal decomposition of nuclear matter |l5l|. In this 
work, by extending the previous work [iq . we study 



nucleon exchange in low-energy nuclear collisions and 
calculate the dispersion of the fragment mass distribu- 
tion 0, 13, 13, S US] ■ Diffusion coefficients for nucleon 
exchange as well as for momentum transfer extracted 
from the SMF approach have the same structure as the 
result familiar from the phenomenological nucleon ex- 
change model. 

The SMF approach is based on a very appealing 
stochastic model proposed for describing deep- inelastic 
heavy-ion collisions and sub-barrier fusion [li, |ll, ]23| . 
In that model, dynamics of relative motion is coupled 
to collective surface modes of colliding ions and treated 
in a classical framework. The initial quantum zero-point 
fiuctuations are incorporated into the calculations in a 
stochastic manner by generating an ensemble of events 
according to the initial distribution of collective modes. 
In the mean-field evolution, couplings of relative motion 
with all other collective and non-collective modes are 
automatically taken into account. In the stochastic ex- 
tension of the mean-field approach, the zero-point (and 
thermal) fiuctuations of the initial state are taken into 
account in a stochastic manner, which is similar to the 
spirit presented in Refs. [H, [3[23|- The initial fiuctu- 
ations are simulated by considering an ensemble of ini- 
tial single-particle density matrices. In this manner, the 
single Slater determinantal description is replaced by a 
superposition of Slater determinants. A member of the 
ensemble, indicated by event label A, can be expressed as 

p\ry,t) = ^<i>*,,(r,t;A)pf^.(ar)<i>,,,(r',i;A), (1) 

ij(7T 

where summations i and j run over a complete set of 
single-particle states <^iar{r,t]\), and a and r denote 
spin and isospin quantum numbers. According to the de- 
scription of the SMF approach [l3|, the element of den- 
sity matrix, p^^ {<Jt) are assumed to be time-independent 

random Gaussian numbers with mean value pfj{<yT) = 
Sijuf^ and the variance of the fluctuating part Sp^AaT) 
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specified by 



-Sjj^Sii'Srr'Saa' [t 



(2) 



Here, nf"^ denotes the average single-particle occupation 
factor. At zero temperature occupation factors are and 
1, and at finite temperature they are determined by the 
Fermi-Dirac distribution. The great advantage of the 
SMF theory is that each Slater determinant A evolves in- 
dependently from each other following the time evolution 
of its single-particle wave-functions in its self-consistent 
mean- field Hamiltonian, denoted by h(p^), according to 
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/l(/)*^.r(r•,^; A). 



(3) 



Following Refs. [21|,l22|, we project the mean-field evo- 
lution on a collective degree of freedom associated with 
nucleon transfer. For the projection, it is useful to intro- 
duce the Wigner distribution for the event A defined as 
a partial Fourier transform of the density matrix. 



(4) 



In this work, we focus on the particular case of head-on 
collisions along the x-axis. We indicate the position of 
the separation plane between the two collision partners 
at X — xq. Then, the mass number of the projectile- like 
fragment in the event A is defined by 
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xo)f\r,p,t). 



(5) 



Other macroscopic variables such as the separation dis- 
tance between fragments and the associated momentum 
can be defined in a similar manner (see [2^1. In the dif- 
fusion model, time evolution of the mass number of the 
projectile-like fragment Ap is described by a Langevin 
equation [25i] . 



dt 



A^P = v{A^P,t)+eAit), 



(6) 



where v{Ap,t) denotes the drift coefficient for nucleon 
transfer. Ignoring memory effects, we consider the quan- 
tity (t) as a Gaussian white noise, which is determined 

with zero mean value Ca(0 = ^ ^^'^ ^ correlation func- 
tion. 



(7) 



where D^a is the diffusion coefficient associated with 
nucleon exchange. In order to extract the diffusion coef- 
ficient, we calculate the rate of change of Ap employing 



the SMF equations. The rate of change of Ap involves 
only the kinetic part of the mean-field Hamiltonian and 
it can be expressed in terms of the reduced Wigner dis- 
tribution on the window as 



dt^^ 



dpx Px J.X, 



p--/ (8) 



where the reduced Wigner distribution f^{x,px,t) is ob- 
tained by integrating over the phase-space variables y, z, 
Py, and pz according to 



f\x,px,t) = 



(9) 



Small fluctuations of the mass number are connected to 
small amplitude fluctuations of the Wigner distribution 
according to 



dt + 



dpx Px c rX/ , \ I 

TT^—Sf {X,px,t)\ 



^Ait). (10) 



In Ref . [16| , we derived an expression for the correlation 
function of the reduced Wigner distribution in the semi- 
classical approximation. Employing the result derived in 
that reference, we have the following expression for the 
nucleon diffusion coefflcient, 



DAAit) 



dpx \Px \ 1 



2TTh m 2 

X >J{/r(a^o,P.,i) [1 - f^^ixo,Px,t)/n(xo,t)] 



(TT 



+ f^^{xQ,Px,t)[l- rp'{xo,Px,t)]/^l{x^,t)]. (11) 

Here r2(a;o, t) is the phase-space volume over the window 
and 



fp'^^rp{xo,Px,t) = JJ dydz J dsr^cxp (^~^PxSx^ 



z, t 



leP/T 



(12) 



is the averaged value of the reduced Wigner distribution 
associated with single-particle wave functions originating 
from the projectile/target. Details on the determination 
of n{xo,t) can be found in Ref. [iGj . 

We note that the expression of the diffusion coeffi- 
cient has the same form as given by the phenomenolog- 
ical nucleon exchange model in Ref. [1^. We also note 
that diffusion coefficients not only for nucleon exchange 
but also associated with other macroscopic variables are 
evaluated in terms of the average evolution specified by 
the standard TDHF evolution. In computations, to em- 
ploy fully quantum mechanical expression for the reduced 
Wigner distribution does not provide a consistent de- 
scription since the diffusion coefficient is derived in the 
semi-classical approximation. A semi-classical form of 
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FIG. 1: (Color online) Evolution of diffusion coefficient calcu- 
lated in SMF approach for ^"Ca+''°Ca (top), ^'^Ni+^'^Ni (mid- 
dle) and ""Zr-l-^'^Zr (bottom) collisions at different center-of- 
mass energies. 

the reduced Wigner distribution can be obtained by ap- 
proximating the Sx dependence of the integrand of the 
expression ([5]) by a Gaussian. The mean value and the 
second moment of this Gaussian are determined by carry- 
ing out a Taylor expansion of the integrand up to second 
order in Sx- 

We carry out calculations for head-on collisions of sym- 
metric '*°Ca-h''°Ca, 56Ni-h56Ni and ^°Zt+^°Zt systems at 
energies just below the Coulomb barrier. Calculations 
are performed using the three-dimensional TDHF code 
developed by P. Bonche and co-workers with the SLy4d 
Skyrme effective force p3] (technical details are given in 
Ref. [21]). Colliding ions approach each other, exchange 
a number of nucleons, and then re-separate. In symmet- 
ric collisions by TDHF, there is no net nucleon transfer, 
i.e., drift is zero. According to the Langevin equation, 
the variance cr\yi{t) — SAp6Ap of fragment mass distri- 
bution, neglecting contributions from the drift term, is 
related to the diffusion coefficient according to i2fl| 

alAt) - 2 / DAAis)ds = iVexc(i). (13) 

In this expression iVexc(i) denotes the accumulated total 
number of exchanged nucleons until time t. The relation 
'^AAi^) — ^cxc{t) follows from nucleon exchange model 
and it was often used to analyze the experimental data. 
In the SMF model, we can calculate the both sides of this 
relation independently. 

Figure [T] illustrates the dependence of diffusion coeffi- 
cients for collision of three different symmetric systems 
at different center-of-mass energies. The Coulomb bar- 



FIG. 2: (Color online) Evolution of 0-^4 calculated in SMF 
approach for ""Ca-t-^^Ca (top), ^'^Ni+^'^Ni (middle) and 
^^Zr-l-^^Zr (bottom) collisions at different center-of-mass en- 
ergies. Number of exchanged nucleons is superimposed by 
the fiUed-circles, fiUed-squares, and fiUed-triangles from high 
to low energies. 



rier energies, which are obtained in the frozen density ap- 
proximation [131, are 54.7 MeV, 103 MeV, and 184 MeV 
for the "^''Ca, ^^Ni, and ^°Zr systems, respectively. The 
magnitude of diffusion coefficient essentially depends on 
the size of the window, the larger the window the larger 
the rate of change of nucleon exchange [1^ . At a given 
center-of-mass energy, diffusion coefficient becomes max- 
imum at the turning point where the size of window is 
the largest. Also, as seen from the figure, due to increas- 
ing overlap of the projectile and target, the magnitude of 
the diffusion coefficient increases with energy. 

Figure [5] illustrates the variances of the fragment mass 
distributions as a function of time for the same sym- 
metric systems at the same center-of-mass energies as 
those in Fig. [TJ Lines are the results obtained by in- 
tegration of diffusion coefficient in Eq. (fTH]) . In each 
case, the corresponding evolution of the number of ex- 
changed nucleons is superimposed by the filled-circles, 
filled-squares, and filled-triangies from high to low ener- 
gies. We also calculate the variances of fragment mass 
distribution using standard TDHF approach. In Table [H 
for the same systems at the same energies, the asymp- 
totic values cr^]3HF(+°°) variances obtained from 
TDHF is compared with the asymptotic values a'^^{+oo) 
calculated in the SMF approach. In the same table, the 
asymptotic values of the number of exchanged nucleons 
Ncxc{+oo) are also given. 

The semi-empirical relation (J^aW — ^exc{t) has 
been extensively used to analyze experimental data 
[28l . [29| . As seen in Table HI the mass variance estimated 
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TABLE I: Asymptotic values of the fragment mass vari- 
ances for *"Ca+'"'Ca, ^^Ni+^^Ni, *Zr+^°Zr collisions in SMF 
((T^^) and TDHF (ctxdhf)- Asymptotic values of the number 
of exchanged nucleons are also given in the last column. 





R fMeVl 

J^c.m. [iv±c V J 


/T^ -i~-.r-r-w^ ( — 1~ 

'-'TDHF V^^y' 




) ' CXC \ ) 


*"Ca+^"Ca 


51.0 


0.004 


0.730 


0.432 




52.5 


0.008 


1.718 


1.441 




53.0 


0.008 


3.790 


3.634 


^'^Ni+^'^Ni 


98 


0.024 


1.288 


0.667 




99.5 


0.088 


2.594 


2.111 




100 


0.142 


5.606 


5.567 


«"Zr+«"Zr 


178 


0.774 


12.98 


14.19 




179 


1.046 


24.09 


27.51 




180 


1.085 


40.40 


41.56 



from SMF is consistent with this relation. Looking at 
the asymptotic values, in all cases, the TDHF results 
cr^j-,jjp(-|-oo) for the variance are much smaller than the 
results cr^^(-l-oo) obtained in the SMF approach and also 
-^exc(+oo)- The failure of TDHF on the description of 
variances of the fragment mass distribution has been rec- 
ognized for a long time as a major limitation of the mean- 
field theory. It appears that the SMF approach cures 
this shortcoming of the mean-field theory. As seen from 
Fig. [2l not only the asymptotic value of <j\j^ but also 
the entire time evolution is very close to the evolution 
of iVcxc (i) • The small differences at energies much below 



the barrier may arise from the estimation of the phase- 
space volume on the window and/or from the Gaussian 
approximation for the reduced Wigner distribution. 

In summary, we investigate the variances of fragment 
mass distributions in heavy-ion collisions at energies near 
the Coulomb barrier employing the microscopic SMF ap- 
proach. By projecting the SMF equation on the mass- 
asymmetry macroscopic variable, we deduce the diffusion 
coefficient associated with nucleon exchange. The expres- 
sion of the diffusion coefficient has a similar structure 
with those familiar from the phenomenological nucleon 
exchange model. Comparison between the calculated 
variance and the number of exchanged nucleon supports 
a strong confirmation for the fact that the SMF approach 
provides a realistic description of dissipation and fluctu- 
ation dynamics at low energies. The stochastic extension 
of the mean-field theory provides a practical solution to 
the estimate of fluctuations of observables at low ener- 
gies. It can be applied not only to nuclear dynamics but 
also for the description of fluctuating dynamics of many- 
body problems in other areas of physics. 
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